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Abstract 



A spin- 1/2 antiferromagnetic(AFM) generalised ladder model is constructed which 
consists of four-spin plaquettes, coupled through weaker exchange interactions, to 
two-spin rungs. In an extended parameter regime, the exact ground state of the 
ladder is determined. In this state, the four-spin plaquettes and the rungs are in 
their ground state spin configurations. In the presence of an external magnetic field, 
the magnetization/site has a plateau structure as a function of the magnetic field. 
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^ ■ 1 Introduction 

O 



Recently, low-dimensional spin systems, particularly the Heisenberg spin lad- 



ders, have been the focus of several analytic, numerical and experimental stud- 
ies [1,2]. Spin ladders can be considered as bridges between one-dimensional 
(ID) and two-dimensional(2D) systems. The ID systems are more or less well 
understood whereas considerable gaps still exist in our understanding of 2D 
systems. The study of ladder models is expected to provide insight on how 
electronic and magnetic properties are modified from ID to 2D. A large num- 
ber of magnetic compounds with ladder structure have been discovered [1,2] 
which exhibit a variety of novel phenomena in the undoped as well as the doped 
states. Different types of spin-1/2 ladder models have been proposed including 
frustrated ladder models and models with modulated exchange interactions [3- 
14]. Ladder models have been studied in zero as well as finite magnetic fields. 
In this paper, we construct a spin-1/2 ladder model with modulated exchange 
interactions. The model consists of four-spin plaquettes connected to two-spin 
rungs. The dominant exchange interactions are within the plaquettes and the 
rungs. The coupling between the plaquettes and the rungs are through weaker 
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Fig. 1. Two-chain ladder model consisting of four-spin plaquettes coupled to 
two-spin rungs. The exchange interaction strengths are as shown in the Figure. 

exchange interactions. Molecular magnets provide another example of spin 
systems consisting of weakly coupled spin clusters[15]. In the case of our lad- 
der model, the spin clusters are the four-spin plaquettes and the two spin 
rungs. 

We show that in a certain parameter regime, the exact ground and low-lying 
excited states of the full ladder model are of the product form, i.e., can be 
written in terms of the exact ground states of the four-spin plaquettes and 
the two-spin rungs. Also, in the presence of an external magnetic field, the 
magnetization/site exhibits the phenomenon of magnetization plateaux. The 
condition for the appearance of a plateau is given by [16] 

S u — m u = integer (1) 



where S u and m u are the total spin and magnetization in unit period of the 
ground state. 



2 Ground and excited states of the ladder model 



The ladder model constructed by us is shown in Fig. 1. The solid lines represent 
the dominant exchange interactions within the four-spin plaquettes and along 
the two-spin rungs. The dotted lines describe the weaker exchange couplings 
between the plaquettes and the rungs. Within a plaquette, the horizontal, 
vertical and diagonal exchange interactions are of strengths J4, J\ and J 3 
respectively. The exchange interaction along a two-spin rung is of strength 
J'. The ladder model consists of alternating four-spin plaquettes and two- 
spin rungs coupled via horizontal and diagonal exchange interactions (dotted 
lines) of strength J2. Periodic boundary condition is assumed to hold true. 
The ladder model generalises a simpler model studied earlier[17] in which the 
horizontal and vertical exchange interactions in a four-spin plaquette are of 
equal strength. The spin Hamiltonian describing the ladder model is 

"<Slj. D 2j + ~Sli+l-S2i+l) 

i=3j+l,j=0,l,- 
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+ ^3(3 n.S 2i+l + O 2i-0 lj+l)] + J' ^ O u.~S 2i 

i=3j,j=0,l,- 

+ ^2i+i) (2) 

i=3j+2,j=0,l,- 

= i^c + H R + i^ci? 

The spin operator 5^, ( S^ 2i ) is associated with the i-th site of the lower 
(upper) chain of the ladder, the site indices are sequential in a chain. The sub- 
Hamiltonians Hq and Hr describe the four-spin plaquettes and the rungs, 
respectively, whereas Hqr contains the exchange couplings between the pla- 
quettes and the rungs. The total spin of each rung is a conserved quantity due 
to the special structure of the Hamiltonian. 

We now determine the ground state of the ladder model using the method 
of 'divide and conquer' [18]. It is easy to show that the state, in which the 
four-spin plaquettes and the two-spin rungs are in their ground state spin 
configurations, is an exact eigenstate of the full Hamiltonian H(Eq. 2). H is a 
sum of three sub-Hamiltonians He, Hr and Hqr- The Hamiltonian Hq + Hr 
acting on the specified state gives back the same state with the eigenvalue, 
E C r, equal to the sum of the ground state energies of all the plaquettes and 
the rungs in the ladder. The sub-Hamiltonian H C r acting on the same state 
gives zero. Thus, the state is an exact eigenstate of H with the eigenvalue 
Ei = Ecr- In an extended parameter regime, the exact eigenstate also turns 
out to be the exact ground state. The proof is as follows:- 

Let Eg be the exact ground state energy of the full Hamiltonian H. Then E g < 
Ei. Let \ip g ) be the exact ground state wave function. Then from variational 
theory, 

E g = £ H 3 |^> + £ H' 3 > J2(E j0 + E' j0 ) (3) 
3 3 j 

j 



where Hj 's are the plaquette Hamiltonians with the ground state energy 
Ej and H'^s are the six-spin cluster Hamiltonians, each of which contains the 
rung exchange interaction Hamiltonian and the eight exchange couplings (four 
horizontal and four diagonal) which connect the rung to nearest-neighbour 
plaquettes. The ground state energy of H] is E' jo . For J 2 < f , E' jo is the 
ground state energy of the rung Hamiltonian. In the ground state, the rung is 
in a singlet configuration. We can now write down the inequality, 

Y,(E io + E'J<E g <Ei (4) 

i 
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Ei is, however, exactly equal to ^2i(E io + E' io ) = Eqr- Thus, E g = Ei, i.e., 
the exact eigenstate of the ladder model is also the exact ground state. The 
exact ground state energy is given by E g = N(E i0 — 3-j) where N is the total 
number of plaquettes as well as rungs in the ladder. 

TABLE I 
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X — yj + J| + J| — Jl J4 — Jl J3 — J3J4 

Table I: The energy eigenvalues and eigenvectors of a four-spin plaquettes with 
exchange interactions of strengths Ji (vertical), J4 (horizontal) and J3 (diagonal). 
^ and ^ are functions of Ji, J 4 and J 3 . 

The eigenvalues and the eigenstates of the four-spin plaquette Hamiltonian are 
shown in Table 1. For the most generalised case, the ground state energy E jo 
is ~( J i+ 2 ? 4+J3) _ x The ground state wave function is of the RVB (resonating 
valence bond)-type. It is a linear combination of two valance bond(VB) states 
with coefficients Ci and c<i depending on the exchange interaction strengths. 
Note that the other singlet-state is also of the RVB-type with the coefficients 
di, d 2 depending on the exchange interaction strengths. Some of the special 
cases of interest are 



4 



(i) Ji = J 4 , J 3 = o 

The ground state energy Ej 

_ — _i 

rf 2 



energy zero and 4 1 = — 1. 



(ii) Ji = J 4 , J 3 + 
For J 3 < Ji, £ i0 = 

301 

(iii) Ji = J 4 = ^3 



.9 7, -I- -h- £i 

ZJ 1 ^ 2 ' c 2 



-2Ji and — = 1. The other singlet state has 

C2 ° 



1. For J 3 > Ji, E'jo = i.e., the other 



RVB state becomes the ground state with ^ = — 1. 



The ground state becomes doubly degenerate. The two states have a pair of 
VBs along either the horizontal or the vertical bonds, with Ejo = — 
(iv) Ji 



The ground state has a pair of singlets along the horizontal bonds with Ej 

3J 4 

(v) 2 J 4 



The ground state has a pair of singlets along the vertical bonds with Ej = 

3Ji 
2 • 

In the generalised as well as the special cases and for J 2 < ^j, the exact ground 
state of the full ladder model is of the product form. The plaquettes are in 
their ground state spin configurations and the rungs are occupied by singlets. 




Fig. 2. Phase diagram of the ladder model (fig.l) in the parameter space of ^ and 



7T for f 



0.75 and J 



0.5. The parameter space below the solid line corresponds 
to the phase in which the exact ground state is a product over the ground states of 
the rungs and the plaquettes. 
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We next checked whether the exact ground state retains its product form 
when J 2 is made larger than For this, we write the total Hamiltonian H 
(Eq. (2)) as a sum over six-spin sub-Hamiltonians, h^s, i.e., H = J2i hi- Each 
sub-Hamiltonian describes a plaquette coupled to a rung. The six-spin sub- 
Hamiltonian can be diagonalised exactly to obtain the ground state energy. 
Again, we use the method of 'divide and conquer'. When the six-spin sub- 
Hamiltonians are added together to obtain the full Hamiltonian, the Ji, J3, J' 
bonds are counted twice and the J 2 bonds only once. We identify the region of 
parameter space in which the full ladder ground state has the product form. 
Fig. 2 shows the phase boundaries, in the parameter space of jj± and ^- for 
fixed values of ^(=0.5) and ^(=0.75). In the parameter regime below the 
phase boundary, the exact ground state has the product form. One finds that 
even for J 2 > 4-, the exact ground state has the product structure. 



3 Magnetization plateaux 



We next include an external magnetic field term — hY^=i SI in the Hamilto- 
nian H (Eq. (2)), where 6N is the total number of sites in the ladder. Let us 
first consider the case of a single 4-spin plaquette in a magnetic field. The mag- 
netic field couples to the z-component of the total spin of the plaquette, S^ ot , 



1/2- 

E 

1/3- 
1/6- 



h h h 

h 

Fig. 3. Plot of magnetization/site m versus external magnetic field h for the 
two-chain ladder model shown in fig. 1. The plot is obtained in the parameter 
region in which the exact ground states in different Sf ot subspaces have the product 
form. Two non-trivial magnetization plateaux occur at m = jr and m = ^. 
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Fig. 4. Phase diagram of the ladder model (fig. 1) in a finite magnetic field and in 
the parameter space of jj- and ^ with ^ = 0.1 and ^ = 0. The regions A, B and 
C are explained in the text. 

which is a conserved quantity. The ground state energy E g (S^ ot ) at h = for 
Sf ot = 0, 1 and 2 can be obtained from Table I. When the external field h ^ 0, 
the ground state in each S^ ot subspace is E g (S^ ot , h) = E g (Sf ot , 0) — hSf ot . The 
ground state magnetization curve can be easily obtained. The magnetization 
per site m is zero from h = upto a critical field h ci . Below h = h Cl , S^ ot of the 
ground state is zero. For h Cl < h < h C2 , S^ ot of the ground state is 1, so that 
m = j and beyond h = h C2 , S^ ot of the ground state is 2, i.e., the saturation 
magnetization m — \ is obtained. Let us first consider the case J 3 < J\. For 
J 3 < J 4 , the ground state in the S^ ot = 1 subspace is ^7) and the critical 
fields are h ci = J 3 + X and h C2 = J x + J4. For J 3 > J 4 , the ground state in 
the Sl ot = 1 subspace is \if) 8 ) and h Cl = J 4 + X and h C2 = J\ + J 3 . At J 3 = J 4 
these two states are degenerate. For J 3 > J 1; as long as J\ < J 4 , the ground 
state in the S^ ot = 1 subspace is with h ci = J\ + X and /i C2 = J 4 + J 3 . For 
Ji > J 4 , the ground state in the S^ ot = 1 subspace is \tps) with h Cl = J 4 + X 
and /i C2 = J\ + J 3 . At Ji = J 4 , |-0 8 ) and j^g) are degenerate ground states. 
Hence depending on the exchange interaction strengths, we get different values 
for h ci and h C2 . 

For the external field h = 0, we have already seen that there is an extended 
parameter regime in which the exact ground state of the full ladder is a product 
of the ground states of the rungs and the plaquettes. We now investigate 
whether the same holds true in the presence of a finite magnetic field. Again, 
we use the method of 'divide and conquer' and the six-spin sub-Hamiltonian 
consists of a plaquette coupled to a rung. For the full ladder, one can identify 




Fig. 5. Phase diagram of the ladder model (fig. 1) in a finite magnetic field and in 
the parameter space of jj- and ^ with ^ = 0.1 and ^ = 0.75. The regions A, B 
and C are explained in the text. 

a region (region A in Fig. 4 and Fig. 5) in parameter space in which for 
< h < h Cl , m is zero. At h Cl , there is a jump in the value of m to m — | 
and a plateau is obtained for h upto h C2 (Fig. 3). When h Cl < h < h C2 , the 
exact ground state has the plaquettes in their S z = 1 ground states and the 
rungs in singlet spin configurations. Since, the number of plaquette is N and 
the total number of sites is 6N, the magnetization/site m in the ground state 
is |. The quantization condition in Eq. (1) is obeyed as unit period of the 
ground state contains six spins so that S u — 3 and the magnetization m u in 
the unit period is 1. At h C2 , the second jump in m from | to |is obtained. 
When h C2 < h < h C3 , the exact ground state has the plaquettes in their S z = 2 
ground states and the rungs in singlet spin configurations. In this case, S u and 
m u in Eq. (1) are 3 and 2 respectively. At h = h C3 , there is a jump in m from | 
to the saturation magnetization |. The value of h C3 , the critical field for which 
the full ladder reaches saturation magnetization is J' + J 2 . There are other 
parameter regions (regions B and C in Fig. 4 and Fig. 5) in the parameter 
space in which the full plateau structure in the m versus h plot, as shown in 
Fig. 3, is not obtained. Fig. 4 shows the phase diagram for the full ladder in 
a magnetic field in the ^- vs. ^ parameter space when J 3 = and ^ = 0.1. 
The region A exhibits the full plateau structure in m vs. h as shown in Fig. 
3. In region B, the jump in m from to | occurs at h — h Cl (Fig. 3) but 
beyond h C2 , the ground state is no longer of the product form. In region C, 
the ground state loses its simple product structure beyond h — h ci . A similar 
phase diagram is shown in Fig. 5, for ^ = 0.75, and with all other parameters 
the same as in the case of Fig. 4. We observe a kink at the point k, where 



8 



J 4 = J3 and a transition from the phase B to phase C occurs. At this point, 
the ground state in S^ ot = 1 subspace changes from \ip s ) to l^). The kink is 
indicative of phase reentrance. When ^- is in the range l.llj^-jl.19, one gets 

the phases B-C-B-C as £ is varied > 0.2). Similarly, for 1.06j J <1.11, 
the phases C-B-C are obtained. 



4 Concluding remarks 

In this paper, we have introduced a generalised ladder model with modulated 
exchange interactions. The model consists of four-spin plaquettes coupled to 
two-spin rungs. In a wide parameter regime, the exact ground state of the 
model is a product over the ground states of the individual plaquettes and 
rungs. In the presence of an external magnetic field, magnetization plateaux 
are obtained when magnetization/site m is plotted as a function of the external 
field h. In an extended parameter regime, the exact ground states in the differ- 
ent magnetization subspaces are of the product form. The generalised model 
includes several other ladder models as special cases. For J2 = J3 = J a = J 
and J' = Ji, the model reduces to the frustrated ladder model introduced 
by Bose and Gayen[19]. As Xian[20] has shown, for — > (y) c ~ 1.40148, 
the exact ground state consists of singlets along the rungs of the ladder. At 
y = (— ) c , there is a transition from the rung dimer state to the Haldane phase 
of the S=l chain. Similar arguments show that for J4 = J3 in our generalised 
ladder model, the ground state in a certain parameter regime is that of a spin 
one chain with modulated exchange interactions. The sequence of exchange 
couplings along the chain has the structure J 4 — J 2 _ J 2 — J a — J 2 — J 2 • 
Other examples of ladder models which are special cases of the generalised 
model, have been given in Section 2. Further studies are needed to obtain the 
phase diagram of the generalised model in the full parameter space. 
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